IHEP 2004-12 



High statistic measurement of the K —> 7r°e v decay 

form-factors 



O.P. Yushchenko, S.A. Akimenko, G.I. Britvich, K.V.Datsko, A. P. Filin, 
A.V. Inyakin, A.S. Konstantinov, V.F. Konstantinov, I.Y. Korolkov, 
V.A. Khmelnikov, V.M. Leontiev, V.P. Novikov, V.F. Obraztsov, V.A. Polyakov, 
V.I. Romanovsky, V.M. Ronjin, V.I. Shelikhov, N.E. Smirnov 1 , O.G. Tchikilev, 

V.A.Uvarov. 

Institute for High Energy Physics, Protvino, Russia 

V.N. Bolotov, S.V. Laptev, A.Yu. Polyarush. 

Institute for Nuclear Research, Moscow, Russia 



Abstract 

The decay K~ — ► 7r°e _ f is studied using in-flight decays detected with the "ISTRA+" 
spectrometer. About 920K events are collected for the analysis. The A+ slope parameter of 
the decay form- factor f+(t) in the linear approximation (average slope) is measured: Alf n = 
0.02774±0.00047(stat) ±0.00032(syst). The quadratic contribution to the form-factor was 
estimated to be A+ = 0.00084 ±0.00027(stat) ±0.00031 (syst). The linear slope, which has 
a meaning of df+(t) / dt\ t=0 for this fit, is A+ = 0.02324±0.00152(stat) ±0.00032(syst). The 
limits on possible tensor and scalar couplings are derived: /t//+(0) = — 0.012±0.021(stat) 
±0.011(syst), fs/f+(P) = -0.0037lo™ (stat) ±0.0041 (syst). 
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1 Introduction 



The decay K — > eun (K e3 ) provides unique information about the dynamics of the strong 
interactions. It has been a testing ground for such theories as current algebra, PCAC, Chiral 
Perturbation Theory (ChPT). The study of this decay has a particular interest in view of new 
two-loop order (0(p 6 )) calculations for decays in ChPT [1, 2]. 

The high-order ChPT calculations make a definite prediction for the quadratic term in the 
vector (/+(£)) form-factor and link the scalar (fo(i)) form-factor linear and quadratic slopes to 
the /+(0) corrections. In turn, /+(0) is known to be crucial for the \ V US \ measurements. In fact, 
the latest measurements do not report any visible non-linearity in the form-factors [3, 4, 5, 6]. 

In this paper we present a high-statistics measurement (~ 919K events) of the Dalitz plot 
density in this decay. The description of the experimental setup, trigger and reconstruction 
procedure can be found in our previous paper [5]. 

2 Selection procedure 

The current analysis is based on the high-statistics data collected during run in Winter 2001. In 
total, 332M events were logged on tapes. This statistics is complemented by about 160M MC 
events generated with Geant3 [7] Monte Carlo program. The MC generation includes a realistic 
description of the setup with decay volume entrance windows, tracking chambers windows, 
chambers gas mixtures, sense wires and cathode structures, Cerenkov counters mirrors and 
gas, the shower generation in EM calorimeters, etc. 

The events with one charged track identified as electron and two or three additional showers 
in the electromagnetic calorimeter are selected for further processing. 

It was observed, that the main background contribution (about 95% of the background 
events) is related to the decay K~ — > 7r~7r°, when the hadronic interaction of the charged 
pion simulates the electromagnetic shower in the calorimeter. Following the method of angular 
selection used in our analysis of the decay [6], we choose the angle between the beam 
particle direction and the vector sum of momenta of the final state track and photons as the 
variable to perform signal-background separation. 

The expected distribution over this angle is shown in Fig.l together with the real data. One 
can observe a clear background peak at small values and a good agreement of the data and 
Monte-Carlo. The cut value was selected in the region where the maximum of the function 
Effsignai/Eff B ack. is reached. 

After that, the momenta of the final state particles were refined with 2C K — > eun kine- 
matic fit. Only the convergence of the fit was required. The missing energy E u = E K — E^ — E v o 
after 2C fit is shown in Fig.2. This variable should exhibit a strong peak at E v ~ in the 
presence of noticeable background contributions. The complete absence of any enhancement 
at small values of the missing energy proofs the quality of the selection procedure. 
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Figure 1: Angle between the beam track 
and the vector sum of final state particles 
momenta. The points with errors are data 
and the solid histogram is MC. The shaded 

area shows the background contribution. 



Figure 2: Energy of the neutrino compared 
with MC. The shaded shows the 
surviving background contribution. 



The signal Monte-Carlo events for Figures 1 and 2 are weighted with the K eS matrix element 
where we use A + = 0.0286 (fixed from our preliminary measurements [5]). 

Finally, 919K events were selected. We estimate the surviving background contribution to 
be around 2.1%. 



3 Analysis 

The most general Lorentz- invariant form of the matrix element for the K~ — > 1~V7T° decay is 
[8]: 

M = ^f^fc^Xl+T 5 )^^ (1) 
2 tuk 

It consists of scalar, vector, and tensor terms. The f± form-factors are the functions of t — 
(Pk — Ptt) 2 - In the Standard Model (SM), the W-boson exchange leads to the pure vector term. 
The scalar and/or tensor terms which are "induced" by EW radiative corrections are negligibly 
small, i.e nonzero scalar or tensor form-factors would indicate the physics beyond the SM. 

The term in the vector part, proportional to /_, is reduced (using the Dirac equation) to 
the scalar form-factor. In the same way, the tensor term is reduced to a mixture of the scalar 
and vector form-factors. The redefined vector (V) and scalar (S) terms, and the corresponding 
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Dalitz plot density in the kaon rest frame (p(E n ,Ei)) are [9]: 

p(E v ,Ei) ~ A ■ \V\ 2 + B • Re(V*S) + C ■ \S\ 2 (2) 
V = f+ + (mi/m K )f T 

S = f s + („H/2m K )f-+(l + 7 ^-^--^-)f T 

\ 2m z K m K m K J 

A = rn K {2E l E v -rn K AE«)-m 2 l {E v --AE«) 
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B = mim K (2E u - AE W ); E v = m K - E t - E n 

yj-> 2 777^ —I— 777^ 

C = m 2 K AE n ; AE n = E™ ax - E n ; E" iax - K 1 n 



2m 



K 



With the selected number of events we can not neglect the V — S interference term proportional 
to the electron mass. The term proportional to /_ is neglected in our analysis. 
For further analysis we assume a general quadratic dependence of /+ on t: 

/+(*) = / + (0) (l + \ + t/m 2 w + \' + t 2 /mi) . (3) 

The procedure of the extraction of the form-factor parameters starts with the subdivision 
of the Dalitz plot region y = 0.12 ^ 0.92; z = 0.55 -r- 1.075 (y = 2E e /m K , z = 2E 7T /m K ) into 
100 x 100 bins. 

The signal MC was generated with the constant matrix element and the amplitude-induced 
weights should be calculated during the fit procedure. One can observe that the Dalitz-plot 
density function p(y, z) in (2) can be presented in the factorisable form: 

P(V, z) = F «( A +> A '+ ? fs, h) ■ K a (y, z), (4) 

a 

where F a are simple bilinear functions of the form-factor parameters and K a (y,z) are the 
kinematic functions which are calculated from the MC-truth information. For each a, the sums 
of K a (y,z) over events are accumulated in the Dalitz plot bins (i,j) to which the MC events 
fall after the reconstruction. Finally, every bin in the Dalitz plot gets weights W a (i,j) and the 
density function r(i,j) which enters into the fitting procedure is constructed: 

r(t,j) = £ F a (\ + , \' + , f s , f T ) ■ W a (i,j) (5) 

a 

This method allows one to avoid the systematic errors due to the "migration" of the events 
over the Dalitz plot due to the finite experimental resolution and automatically takes into 
account the efficiency of the reconstruction and selection procedures. 

To take into account the finite number of MC events in the particular bin and strong 
variation of the real data events over the Dalitz plot, we minimize a — C function defined as 
[10]: 



-C = 2J2 rij In 



rj V rrij + l t 



2 5> 3 -+m 3 - + l)ln| - . g | , ( 6 ) 



1 + u. 

m 3 
n i 

rrij+l 



3 



where the sum runs over all populated bins, and rij, rj and rrij are the number of data events, 
expected events and generated Monte Carlo events respectively. For large rrij Eq. (8) reduces 
to the more familiar expression 

-C = YH 2 ( r i ~ %) + 2n j lnn i/ r i] 

3 

The radiative corrections were taken into account by re- weighting every Monte-Carlo event, 
using MC-truth information, according to the recent calculations in [11]. 

The minimization is performed by means of the "MINUIT" program [12]. The errors are 
calculated by "MINOS" procedure of "MINUIT" at the level AC = 1, corresponding to 68% 
coverage probability for 1 parameter. 

4 Results 

A fit of the K e3 data with f s = f T = \' + = gives A+ = 0.02774 ± 0.00047. The total number 
of bins is 6991 and x 2 / n df = 0.97. The quality of the fit is illustrated in figures 3 and 4 where 
the projected variables y = 2E e /m K and z = 2E 7T o/m K are presented. 




Figure 3: Y distribution. Figure 4: Z distribution. 

The points with errors are the real data The points with errors are the real data 
and the shaded area - signal MC. and the shaded area - signal MC. 

The Table 1 represents fits with possible nonlinear term in / + (Eq. 3) as well as the fits 
with tensor and scalar contributions (Eq. 1). 

Every row of the Table 1 represents a particular fit where the parameters shown without er- 
rors are fixed. To qualify the statistical significance of the particular fit, we also show the change 
in the \ 2 value obtained with respect to the fit without non-linear or anomalous contributions. 
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The second row shows a fit where the nonlinearity is allowed in f+(t). One can observe 
A + — A^ (Fig. 5) correlation that results in a significant A + errors enhancement and visible shift 
of A + parameter. 
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0. 


0. 
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0. 


-0.012 ±0.021 


0. 


-0.3 


0.02771 ± 0.00047 




0. 


0. 


-0.0059l°S 


-0.5 


0.02324 ±0.00152 


0.000^ 


54 ±0.00027 


-0.012 ±0.021 


0. 


-9.9 


0.02325 ±0.00152 


0.000^ 


54 ±0.00027 


0. 


U - UU <J' -0.0056 


-9.9 



Table 1. The K e ^ fits. 



x 10 
0.14 p 
0.13 r 




0.04 h 1 

0.02 0.021 0.022 0.023 0.024 0.025 0.026 



Figure 5: la and 2a contours in the A + — A+ plane. 

. We should consider the A + parameter in the linear approximation as a mean slope over all 
the physical g 2 -region, while the true value of the linear slope (df+(t) / dt\ t =o) should be taken 
from the non-linear fit. 
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We also perform model-independent one-dimensional (y) fits where the data in every of the 
100 q 2 /m 2 bins were fitted independently. The resulting distribution is shown in Fig.6. The 
normalization /+(0) = 1 is assumed. The visible non-linearity can be observed in Fig. 7, where 
the ratio / + (i)// + (0)/(l + \ + q 2 /m%) is presented. The parabolic curve represents the fit with 
the quadratic non-linearity in the form-factor. 
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Figure 6: The value of f + (t)/ f+(0) obtained 
in the model-independent fits. 



Figure 7: The value of 
/+(f)//+(0)/(l + A + g 2 /0- The fit with 
non-linear contribution is shown. 



This non-linearity can not be explained by a possible scalar contribution (that also results 
in the enhancement of the number of events at large values of q 2 ). The row 4 of the Table 
1 represents a search for the scalar term with the vector form-factor set to be linear. The 
resulting value of /s//+(0) is compatible with zero. 

We also perform a model-independent fit to extract simultaneously f+(t) and fs(t)- The 
resulting distribution for the value fs(t)/f+(0) is shown in Fig. 8. The value of the scalar 
contribution is compatible with zero with strong enhancement of the errors at small values of 
t. This enhancement is explained by the dependence of the scalar contributions (Eq. 2) on the 
Dalitz variables. One can observe that the leading term \S\ 2 is proportional to t and vanishes 
at t -> 0. 

The last row of the Table 1 represents a fit with both scalar contribution and the quadratic 
term in the vector form-factor. 

We also do not see any tensor contribution in our data (rows 3 and 5 in the Table 1). 
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Figure 8: The value of fs(t)/f + (0) obtained in the model-independent fit. 

Different sources of systematics are investigated. We allow variations of the electron selection 
cuts and angular cut. The Dalitz plot binning, signal and background MC variations are also 
applied. 

The systematic errors are summarized in the Table 2. We can conclude that the main 
contribution into systematic error comes from the main selection cut variation and from the 
limited amount of the background MC. 



Source 


A+ 


a; 


/t//+(0) 


fs/U(0) 


selection 


0.00017 


0.00013 


0.005 


0.0020 


angular cut 


0.00020 


0.00021 


0.006 


0.0020 


Dalitz plot binning 


0.00004 


0.00006 


0.001 


0.0005 


Signal MC variation 


0.00006 


0.00004 


0.002 


0.0006 


Backg. MC variation 


0.00016 


0.00018 


0.008 


0.0026 


Total 


0.00032 


0.00031 


0.011 


0.0041 



Table 2. The systematic error contributions. 



7 



5 Summary and conclusions 

The K~ 3 decay has been studied using in-flight decays of 25 GeV K~ detected by the "ISTRA+" 
magnetic spectrometer. 

The A + parameter of the vector form-factor in the linear approximation (average slope) is 
measured to be: 

Af = 0.02774 ± 0.00047 (stat) ± 0.00032 (syst). 

A visible non-linear contribution is observed for the first time : 

A+ = 0.00084 ± 0.00027 (stat) ± 0.00031 (syst). 

With the quadratic term in the vector form-factor the linear slope (which has a meaning of 
df + (t)/dt\ t=0 ) is determined to be: 

A + = 0.02324 ± 0.00152 (stat) ± 0.00032 (syst). 

The limits on possible tensor and scalar couplings are derived from the fit with quadratic vector 
form-factor: 

f T /f + (0) = -0.012 ± 0.021 (stat) ± 0.011 (syst); 
fs/f+(0) = -0.00371S (stat) ± 0.0041 (syst). 

The value of the quadratic term in the vector form-factor is in a good agreement (within 
the statistical errors) with the 0(p 6 ) ChPT predictions [1, 2]. 

The obtained limits on the scalar and tensor terms can be used to get limits on several 
exotic models. 

The leptoquark- induced amplitudes for K e3 decay were considered in [13]. In this model the 
relation between tensor and scalar terms is fixed: fr = — 6.87 fs- The fit of our data with this 
constraint gives the following scalar contribution: /s//+(0) = 0.0046to;oogy. This result can be 
converted in the upper limit |/s// + (0)| < 0.011 (90% C.L.). Using the expression in [13]: 

fs = V% m 2 K -ml 1 
f + (0) 16G F \V SU \ (m s -m u )m K Al Q ' 1 ' 

where Alq is the ratio of leptoquark mass to the square of the Yukawa-like coupling, we get 

A LQ > 2.55 TeV (90% C.L.) 
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